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Summary 

Two analytical models, one based on simple hydro- 
dynamic lubrication (HL) and the other on elasto- 
hydrodynamic lubrication for materials of low elastic 
modulus (soft EHL) are presented and compared to 
delineate the dominant physical parameters that govern 
the mechanics of a gaseous film between a small drop of 
lubricant and the outer race of a ball bearing. Both 
models are based on the balance of gravity forces, air 
drag forces, and air film lubrication forces and 
incorporate a drag coefficient C D and a lubrication 
coefficient C L to be determined from experiment. The 
soft EHL model considers the effects of droplet 
deformation and solid surface geometry, whereas the 
simpler HL model assumes that the droplet remains 
essentially spherical. The droplet’s angular position 
depends primarily on the ratio of gas inertia to droplet 
gravity forces and on the gas Reynolds number and 
weakly on the ratio of droplet gravity forces to surface 
tension forces (Bond number) and geometric ratios for 
the soft EHL. 

An experimental configuration in which an oil droplet 
is supported by an air film on the rotating outer race 
of a ball bearing within a pressure-controlled chamber 
produced measurements of droplet angular position as a 
function of outer -race velocity, droplet size and type, and 
chamber pressure. For Bond numbers less than 0.5 
(larger surface tension forces), both models correlate the 
experimental data with a nearly constant C D of 0.4 and a 
C L of 3.0 for HL or 1.0 for soft EHL. For larger 
variations in Bond number B and Reynolds number Re/), 
the soft EHL model shows superior correlation of the 
experimental data with experimentally determined C D 
versus Re/) and C/ versus B. 


Introduction 

Although the amount of lubricant needed to sufficiently 
lubricate a nonconformal contact elastohydrodynamically 
is small, a critical amount is needed to maintain a fully 
flooded conjunction. Nonconformal contacts are 
normally found in mating gear teeth, cams and followers, 
and rolling-element bearings. If the amount of fluid in 
the conjunction goes below the critical amount, the 


conjunction becomes starved and the minimum film 
thickness is considerably reduced, as shown by Hamrock 
and Dowson (1977). If the starvation is continued, the 
nonconformal contact will eventually fail. 

Kingsbury (1969) reported the failure of large lubricant 
droplets to wet and enter the conjunction of a rolling- 
element bearing, thereby causing lubricant starvation. 
Starvation is caused by the presence of a gaseous film 
barrier between the droplet and the bearing race surface. 
The present work presents experimental and theoretical 
results that delineate the physical parameters that govern 
the behavior of a stable droplet on a gas film. 

The experimental apparatus consisted of an outer race 
of a ball bearing rotating in the presence of air within a 
pressure chamber while the equilibrium position of 
various liquid droplets was measured. The droplet 
position was recorded as the angle subtended by the 
radial line from the center of rotation to the droplet 
center and the vertical line through the rotational center, 
for the range 20° to 70°. 

The physical parameters and their ranges were 

(1) Droplet radius, 1 to 1.75 mm 

(2) Velocity of the outer race, 2 to 10 m/s 

(3) Chamber pressures, 98.6, 68.9, and 34.5 kPa 

(4) Three lubricating oils 

Besides accurately measuring the position of the 
droplet for these parameters, an analytical model was 
developed that describes the position of the droplet as a 
function of the forces acting on it. The forces considered 
were gravity, drag, skin friction, and lubrication. Two 
approaches to quantifying the lubrication force were 
considered and compared: a semiempirical hydro- 
dynamic model (HL) that assumes the droplet shape to be 
essentially spherical, and an elastohydrodynamic model 
for materials of low elastic modulus (soft EHL), which 
accounts for droplet deformation. By weighing the 
relative importance of the physical parameters in the 
problem, the effects of surface tension, gravity, air 
density, and air viscosity were delineated through the 
dimensionless ratios of the Reynolds and Bond numbers 
and the ratio of air inertia to oil gravity forces (a 
modified Froude number). The predictions of the droplet 
angular position versus modified Froude, Bond, and 
Reynolds numbers were compared for the two models. 
These results were compared with the experimental 
measurements taken over a wide range of operating 
parameters. 
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Symbols 


B 

Bond number, p 0l fcig/a 

Qd 

drag coefficient 

C EHL 

elastohydrodynamic lubrication coefficient 

C L 

C L , 2 /{C l> i)1/2 

C L, 1 

hydrodynamic lubrication coefficient 

c l, 2 

skin friction coefficient 

D x 

diameter of contact ellipse along direction of 
motion, m 

D y 

diameter of contact ellipse along transverse 
direction, m 

d d 

diameter of droplet, m 

E 

modulus of elasticity, Pa 

E' 

effective modulus of elasticity, 
2/[(l-v 2 d )/E d +(l-p 2 d )/E b ], Pa 

S 

complete elliptical integral of second kind 

F 

modified Froude number, or EHL lubrication 
force, N 

8 

gravitational constant, m/s2 

h 

film thickness, m 

^min 

minimum film thickness, m 

k 

ellipticity parameter, D y /D x 

Pc 

chamber pressure. Pa 

P 0 

ambient pressure, Pa 

R 

curvature sum, m 

R c 

correlation coefficient 

Re/) 

air Reynolds number, p a [ v Vd d /rj aiT 

r 

radius of curvature, m 

r d 

radius of oil droplet, m 

u 

surface velocity in direction of motion, 
(u d + K)/2, m/s 

V 

velocity of outer race, m/s 

We 

Weber number, p^[ r V 2 d d /a 

w 

weight of oil droplet, N 

V 

absolute viscosity at gauge pressure, N s/m 2 

9 

angular position of droplet 

V 

Poisson’s ratio 

P 

lubricant density, N s 2 /m 4 

a 

surface tension, N/m 


Subscripts: 
b outer race 

d oil droplet 

x,y coordinate system 

0 ambient conditions 

Experimental Apparatus 

An apparatus was assembled to exhibit the phenom- 
enon of air film droplet suspension on a rotating surface 
and to make the characteristics of this phenomenon 
amenable to measurement. The equipment (fig. 1) 
consisted of the outer race of a ball bearing mounted with 
its axis of rotation horizontally supported by two ball 
bearings within a pressure chamber. A removable glass 
viewing port was placed normal to the axis of rotation 
and close to the droplet location for ease of droplet 
insertion. The outer race and shaft assembly were driven 
by a variable-speed motor drive through an O-ring seal in 
the chamber opposite the viewing port. The rotational 
speed of the motor was proportional to the input voltage, 
measured by a digital voltmeter. A strobescope was used 
to obtain precise rotational speeds. Suitable couplings 
were placed between the sealed shaft of the motor and the 
outer-race shaft assemblage to ensure vibration-free 
rotation. The three chamber pressures (98.6, 68.9, and 
34.5 kPa) were measured by an absolute pressure gauge 
and were maintained by adjusting a bleed valve to the 
atmosphere while evacuating the chamber with a vacuum 
pump. 

The droplet’s size and angular position on the race 
were determined with a cathetometer mounted on a 
horizontal traversing compound. A sectional view of the 
outer race and a typical droplet location are shown in 
figure 2. 


Test Procedure 

Before each day’s tests the cathetometer was leveled 
and adjusted so that the telescopic site moved in a plane 
parallel to the plane of rotation of the outer race. The 
eyepiece scale of the cathetometer was calibrated by 
placing a steel ball of precisely 0.318 cm (0.125 in) 
diameter on the outer race and measuring its diameter as 
a given number of units of the eyepiece scale. The 
horizontal and vertical readings of the position of the 
outer-race rotational center were checked frequently so 
that they, combined with the droplet’s horizontal and 
vertical position, could be used to determine the droplet’s 
angular position 6. 
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Figure 1.— Experimental apparatus, showing rotating outer race, pressure chamber, drive motor, cathetometer, and traversing compound. 



Figure 2. — Sectional view of outer race and drop location. 


The ambient pressure in the laboratory was compared 
with the absolute pressure gauge’s reading to ensure 
proper operation of the gauge. The ambient temperature 
of the laboratory was checked to ensure that the labora- 
tory’s constant-temperature control was functioning. 

Once these preliminary steps had been followed, the 
test oil was drawn into a small syringe, and a droplet was 
inserted onto the rotating outer race (fig. 3). This was 
accomplished with the window removed from the 
chamber and at ambient conditions. Initially the 
rotational velocity for a stable droplet angular position 


was not known. After a few trials the approximate 
velocity range for droplet insertion was determined. The 
droplet’s size could be altered by adding fluid to, or 
extracting fluid from, a droplet while it was on an air film 
in the apparatus. Once a droplet of approximately the 
desired size was situated on the rotating outer race, the 
glass viewing port was placed on the O-ring and held in 
place by rubber bands. 

The droplet diameter was measured in the tangential 
and radial directions while the droplet was suspended at a 
stable position on the rotating race, usually at about 45°. 
Only for the largest droplets tested was there a 
measurable difference between the tangential and radial 
diameters, since the inaccuracy in the measurements of 
size on the fluctuating droplet was of the same order as 
the variation in droplet dimensions. The average of the 
two readings was used throughout the runs on a given 
droplet, although at angles near 60° the droplet was 
visibly misshapen by the excessive air inertial forces. 

The normal procedure was to work the droplet from 
the lowest angle (found by trial and error to be about 15°) 
through to the highest angle for which the droplet 
appeared stable by increasing the velocity of the race 
systematically at ambient pressure. At about 60° to 65° 
the droplet stood off the race and began to oscillate in the 
axial as well as the tangential direction. The droplet’s 
position was recorded at each rotational velocity, usually 
in increments of 100 rpm. With the same droplet, if the 
experimenter had avoided losing the droplet at either end 
of the angle domain, the pressure was dropped to 
68.9 kPa while the velocity of the race was manually 
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Figure 3. — Droplet insertion onto rotating outer race with hypodermic syringe. 


raised to maintain the droplet in the most stable angle 
range of about 45°. The race velocity was again raised 
from the lowest “safe” value through to the highest value 
while angles were recorded and the pressure was checked 
to ensure constancy. The pressure was then lowered to 
34.5 kPa while the droplet’s angular position was 
preserved at 45° by manually raising the race velocity. 
The reproducibility of this overall procedure for a given 
droplet was better than 1 to 2 percent and from droplet to 
droplet was of the order of 5 percent since the most 
imprecise measurement stems from the droplet size 
determination. 


Model 

A simplified view of a droplet of liquid positioned on 
the rotating outer race of a ball bearing looking along an 
axis through the droplet center and parallel to the axis of 
rotation is sketched in figure 4. The important forces are 
indicated in the figure. The tangential and radial force 
balance can be written for a simple hydrodynamic model 
as 


CpPair ^ + C Lt 2 r?a ir V^h _ w sin e= ! 
2 h 



Figure 4.— Schematic of droplet on outer race with characteristic forces 
of weight, air drag, and air film lubrication. 


where p air and r/ a j r are the air density and viscosity, 
respectively, w is the oil droplet weight, r d is the droplet 
radius, h is the film thickness, V is the outer-race 
velocity, 0 is the droplet angular position, C D is the drag 
coefficient, and C L 2 and C L \ are skin friction and 
lubrication coefficients, respectively, introduced to 
absorb geometric considerations not included in this 
simplified model. The frontal area of the droplet and the 
droplet contact area are both set to 7rr^ with the 
coefficients C D , C L \, and C L2 expected to correct for 
the inadequacies of these assumptions. 
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By combining equations (1) and (2) and rearranging 
terms, the following relation for predicting 0 was 
obtained: 


sin 0 = 


3C ZJPair^ 1 2 | r ( 37/ air F \ 1/2 (COS0) 1/2 

Spoilt Z ' , 2 Mp oil g^/ (Q,.i ) 172 


(3) 


where p 0 ij is the oil droplet’s density and g is gravity. 

In figure 4 the droplet is shown as spherical, implying 
that the air-oil surface tension forces are large as com- 
pared with gravity and air inertia. If in fact the droplet 
shape deviates appreciably from spherical because surface 
tension forces are not dominant, a more complex model 
than that just presented is required to predict 0 as a 
function of fluid properties, race velocity, and gravity. 
The elastohydrodynamic theory of Hamrock and Dowson 
(1978) for materials of low elastic modulus (soft EHL) is 
an attractive candidate for a model that considers 
deformation of the droplet in spite of the fact that 
droplet deformation cannot be precisely characterized as 
deformation of an elastic solid. The lubrication force in 
equation (2) therefore is represented by the elasto- 
hydrodynamic lubrication model. Furthermore in 
equation (1) the skin friction term can be more accurately 
described by defining the contract area of the deformed 
drop as i rD x D y /4 instead of tt as was done in the 
hydrodynamic model. The tangential and radial force 
balances that result for the soft EHL model are 


CdP^V 2 ^ + C EHLVair V*D x D y _ wsin g_ Q (4) 

^ ™min 

£ '^[ 7 - 43 (^) <, - o - 85e ' o - 3 ' t) (ifJ 0 “] 4K 

-w cos 0 = 0 (5) 


where 


r _ r d r b,x 
X r d + r b,x 


Following the convention that convex surfaces have 
positive curvatures and concave surfaces, negative 
curvatures, the droplet will have a positive radius and the 
two radii of the outer race will be negative. 

Surface tension and gravity affect the droplet shape: 
the local radii of curvature of the droplet near the contact 
area must be such that the surface tension forces can 
support the hydrostatic pressure within the droplet. This 
internal pressure is characterized as p 0 \\d d g cos 0, where 
the droplet height is assumed to be nearly d d . The low 
modulus of elasticity of the droplet, approximated as the 
pressure within the droplet resulting from the surface 
tension forces, is then also equal to p 0 [\d d g cos 0. 
Furthermore the modulus of elasticity of the droplet is 
orders of magnitude less than that of the outer race, and 
Poisson’s ratio of the liquid is assumed to be 0.5. Making 
use of this, the effective modulus of elasticity becomes 


, _ 8p oi igrf rf cos 6 

3 


E' = 


From Hamrock and Brewe (1983) simplified equations 
for the ellipticity parameter and the elliptic integral of the 
second kind are 


*-M *r 

-•(*)(¥) 


These equations can be used to define the elliptical 
contact diameters in the x direction (direction of motion) 
and in the y direction (transverse direction), which are 
expressed in equation (4) as 


/ 6Swi? \ 1/3 
V irkE' ) 


2 / 6k 2 &wR ^ 1/3 


/ 


j? ~ b>y 
y r d +r bty 


u — 


u d +V 

2 


V_ 

2 


where 

1 _ 1 1 
~R~~ R~ x + Ry 


1 -*d , \zA 

E d E b 


and 
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The radial force balance can be simplified and put into a 
form similar to that of the hydrodynamic model. 


Vmr Yl!^ \ 

^min ' ' ^min ' 


2.96(1 -0.85 e” 0 * 31 ^ 1 * 54 



- w cos 0 = 0 (6) 


where Fis the modified Froude number and C d (Rq d ) and 
C l (B ) are relations that we obtain by comparing theory 
and experiment for a variety of fluids. Note that the 
product of the first and third groups in equation (9) is 1/2 
Weber number, p a \ T V 2 d d /a, which plays a role in the 
droplet shape if air inertial forces become excessively 
large as compared with the surface tension. Dimen- 
sionless scalings with the ranges experienced in the tests 
are presented in table I. 


The soft EHL model includes droplet and outer -race 
geometry and thus eliminates the need for the empirical 
coefficient C L i of the hydrodynamic model. More 
significantly, the r d /h m in factor magnifying the ratio of 
normal force to tangential force is to the 0.54 power as 
compared with the first power in the hydrodynamic 
model. 

By combining equations (4) and (6) and rearranging 
terms, the following relation for predicting 6 is obtained: 


sin 6 = 


3Q*>air^ 2 

Spoiled 



(1 -0.85 e-0.3U)4^ 


(7) 


TABLE I.— DIMENSIONLESS SCALINGS 


Dimensionless scalings 

Low 

High 

Air inertia p a K 2 / 2 P a V 2 r d 1 

= We^ 

Surface tension 2 a/r d 4 a 8 

0.1 

1.5 

Gravity P 0 a^r d g p oiI gr d n 

0.3 

1.4 

— — — D 

Surface tension a 

Air inertia p air V 1 /! p air Vd d 1 

= = - Re D 

Air viscosity r/ air V/r d 4r/ air 4 

50 

500 

Air inertia p air K 2 /2 p air Vh 1 ( h \ 

Air film viscosity rj air V/h 2 7j air 4 \r d / 

1.5 

40 

Air film inertia P^V 2 /Ar d p^ x Vd d ( h V 
Air film viscosity i? air V/h 2 \r d J 

0.01 

0.5 


Comparing equations (3) and (7) shows the effects of 
the soft EHL model primarily in the power of the 
parameter 3?7 air F/4p oil g/^ being 0.35 instead of 0.5 and in 
the power of the cos 0 term being 0.65 instead of 0.5. 

Equations (3) and (7) are both of the form 


p_pfpair ^ 2 VairV ^ ^ x i^l 


( 8 ) 


Results 

Three oils were tested extensively: synthetic paraffin 
(XRM 109 F-3), fluorosilicone (XF 10301), and four-ring 
polyphenyl ether (MCS 210). The relevant properties of 
these fluids at 25 °C are given in table II. 


where x/r d represents the length ratios r bx /r d and 
r b,y/ r d* The drag coefficient C D is expected to be a weak 
function of the air Reynolds number P a [ r Vd d /r} a [ T , which 
is twice the ratio of the first and second groups in 
equation (8). The lubrication coefficient C L might have 
some dependence on surface tension for large droplets. 
Comparing gravity forces with surface tension results in 
the Bond number Pow^g/a. Equation (8) can be put in the 
alternative form 

e=f^F,Rt D ,B, (9) 


TABLE II.— FLUID PROPERTIES AT 25 °C 


Fluid 

Surface 

tension, 

N/m 

Density, 

kg/m 3 

Viscosity, 
kg/m s 

Synthetic paraffin 
(XRM 109 F-3) 

0.0303 

840 

0.840 

Fluorosilicone 
(XF 10301) 

.0250 

1150 

.055 

Four-ring poly- 
phenyl ether 
(MCS 210) 

.0448 

1180 

.177 
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(a) V=2.85 m/s; We D = 0.47; R ez) = 410. 

(b) V= 5. 30 m/s; We/>=1.6; Re D = 760. 

(c) K= 7.74 m/s; We D = 3.5; Re^llOO. 

Figure 5.— Four-ring polyphenyl ether (MCS 210) droplet of 1.08-mm radius in air at Bond number of 0.30 and chamber pressure of 98.6 kPa 
for varying outer-race velocities and Weber and Reynolds numbers (based on diameter). 


The radii of curvature of the outer race are 
r b, x = "38.9 mm and r by = - 6.55 mm, and the rms 
roughness is approximately 0.5 to 1 /un in the y direction 
and 0.1 to 0.2 /mi in the x direction. 

Figures 5 and 6 illustrate the characteristics of a 
1.08-mm-radius droplet of four-ring polyphenyl ether 
(MCS 210) at a small Bond number (0.30). The velocity 
of the outer race was raised from 2.85 to 5.30 to 7.74 m/s 
at a chamber pressure of 98.6 kPa (fig. 5). The dom- 
inance of the surface tension forces gave the droplet a 
generally spherical appearance. The Weber and Reynolds 
numbers are given on each figure. The maximum We of 
3.5 (fig. 5(c)) occurred at the highest velocity and largest 
angle, near 60°, but 1/8 We = 0.40, which from table I is 
most nearly a measure of air inertia to surface tension 
forces, is still less than unity and the droplet did not 
appear to deviate substantially from spherical. 

Figure 6 shows the effect of reducing the chamber 
pressure from 98.6 to 68.9 to 34.5 kPa at a substantially 
constant velocity near 8 to 9 m/s for the same droplet of 


MCS 210 at a Bond number of 0.30. Since the air density 
in the chamber was proportional to the chamber 
pressure, 6 fell from 60° to 50° to about 40° as the 
pressure was reduced. 

Figures 7 and 8 illustrate the alterations in the shape of 
a 1.98-mm-radius droplet of MCS 210 at a Bond number 
of 1.0. Raising the outer -race velocity from 5.30 to 
7.54 m/s at a chamber pressure of 98.6 kPa (fig. 7) 
caused the Weber number to vary from 3.0 to 6. 1 to 10.0. 
In figure 7(c), 1/8 We exceeds unity and the shape of the 
droplet is considerably different from a sphere. Figure 8 
again shows the effect of chamber pressure reduction for 
this large drop of MCS 210 at a fixed outer -race velocity 
of 9.8 m/s as 0 decreased from 50° to 40° to about 15°. 

The results of 6 versus outer-race velocity V for an oil 
droplet of known physical properties and radius r d at a 
given chamber pressure P c were used in conjunction with 
the prediction equations (3) and (7) to determine the drag 
and lubrication coefficients. Equation (3) was put in the 
form 
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(a) K= 7.74 m/s; We D = 3.5; Re D =1100; P c = 98.6 kPa. 

(b) K=7.74 m/s; We D = 2.4; Re Z) = 780; P c = 68.9 kPa. 

(c) V= 8.96 m/s; We D = 1 .6; Re D = 450; P c = 34.5 kPa. 

Figure 6. — Four-ring polyphenyl ether (MCS 210) droplet of 1.08-mm radius in air at Bond number of 0.30 for varying outer-race velocities, 
Weber and Reynolds numbers (based on diameter), and chamber pressures. 
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(a) V= 5.30 m/s; We D = 3.0; Re o =1400. 

(b) K=7.54 m/s; We fl = 6.1; Re D = 2000. 

(c) K=9.78 m/s; We £) = 10.0; Re D = 2600. 

Figure 7.— Four-ring polyphenyl ether (MCS 210) droplet of 1 .98-mm radius in air at Bond number of 1 .0 and chamber pressure of 98.6 kPa for 
varying outer-race velocities and Weber and Reynolds numbers (based on diameter). 
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(a) We 0 = 10.0; Re D = 2600; P c = 98.6 kPa. 

(b) We o = 7.0; Re c = 1800; P c = 68.9 kPa. 

(c) We£> = 3.5; Re z> = 900; P c = 34.5 kPa. 

Figure 8. Four-ring polyphenyl ether (MCS 210) droplet of 1.98-mm radius in air at Bond number of 1.0 and outer-race velocity of 9.78 m/s 
for varying Weber and Reynolds numbers (based on diameter) and chamber pressures. 
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y— Cipc+Cj) 


( 10 ) 


where 


Spoilt sin 6 
3p air F2 


a (’/air VPo\\8 COS 0/ 3) 1/2 


and 


Pair — 


_ Poil^c 


oil 


which yielded C L , Cp, and the correlation coefficient by 
linear regression. 

A similar approach with equation (7) yielded C L , C D , 
and the correlation coefficient for the soft EHL model. 
The coefficients for the two models for the three fluids 
and three chamber pressures (table III) suggest that the 


soft EHL model is superior to the simple HL model. The 
soft EHL model showed its superiority at droplet sizes 
and angular positions where the droplet shape deviated 
substantially from spherical. When the higher angle data 
were eliminated from consideration in determining the 
drag and lubrication coefficients from the hydrodynamic 
equation (10), higher correlations were obtained. 

Both models’ correlation coefficients improved for low 
Bond and Weber numbers. These conditions are gener- 
ally realized for small droplet radii and low chamber 
pressures (i.e., low air density, table III). 

The drag coefficients for the soft EHL model shown in 
table III and for intermediate radii not included in table 
III as a function of Reynolds number for all three oils 
and pressures are shown in figure 9. The range of 
Reynolds number and uncertainty in the drag coefficient 
is indicated by the bars. The results were below those of 
wind tunnel data on smooth spheres. The scatter in the 
drag coefficient results did not correlate with the 
variation of fluid type. 

The lubrication coefficients for the soft EHL model 
shown in table III and for intermediate radii not included 
in table III are shown graphically as a function of Bond 
number in figure 10 for the three fluids at the three 


TABLE III.— COMPARISON OF HYDRODYNAMIC (HL) AND ELASTOHYDRODYNAMIC (EHL) LUBRICATION 

MODELS FOR THREE FLUIDS AND PRESSURES 


Fluid 

Viscosity 

ratio, 

^oil^air 

Chamber 

pressure, 

Pc 

kPa 

Droplet 

radius, 

r d * 

mm 

Bond 

number, 

B 

Reynolds 

number, 

Re 0 

Lubrication 

coefficient, 

C L 

Drag 

coefficient, 

c D 

Correlation 

coefficient, 

R c 

HL 

EHL 

HL 

EHL 

HL 

EHL 

Synthetic paraffin 

4.600x 10~“ 

98,6 

1.06 

0.300 

460-860 

2.9 

0.91 

0.36 

0.38 

0.987 

0.998 

(XRM 109 F-3) 



1.75 

.820 

1135-1900 

1.8 

.63 

.34 

.34 

.864 

.994 



68.9 

1.06 

.300 

360-725 

3.1 

1.0 

.37 

.38 

.979 

.996 




1.75 

.820 

930-1530 

2.0 

.66 

.36 

.35 

.905 

.995 



34.5 

1.06 

.300 

200-445 

3.1 

.98 

.44 

.47 

.997 

.999 




1.75 

.820 

600-800 

1.7 

.45 

.43 

.45 

.972 

.999 

Four-ring poly- 

0.980x10'“ 

98.6 

1.16 

0.347 

525-1130 

3.5 

1.1 

0.43 

0.43 

0.970 

0.995 

phenyl ether 



1.69 

.735 

1280-2000 

1.6 

.48 

.39 

.39 

.919 

.998 

(MCS210) 














68.9 

1.14 

.336 

430-950 

3.3 

1.0 

.39 

.40 

.982 

.997 




1.69 

.735 

1020-1530 

1.7 

.47 

.41 

.42 

.993 

1.000 



34.5 

1.10 

.313 

180-500 

3.3 

.96 

.43 

.48 

.997 

1.000 




1.69 

.735 

600-770 

2.5 

.67 

.42 

.44 

.991 

.999 

Fluorosilicone 

0.350x10'“ 

98.6 

1.00 

0.450 

490-870 

2.7 

0.80 

0.43 

0.43 

0.992 

0.999 

(XF 10301) 



1.74 

1.367 

1130-1980 

.78 

.33 

.44 

.42 

.611 

.993 



68.9 

1.00 

.450 

380-685 

2.8 

.83 

.45 

.47 

.997 

1.000 




1.74 

1.367 

925-1590 

1.2 

.38 

.44 

.44 

.824 
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Figure 9. — Drag coefficient as function of Reynolds number (based on 
diameter) for four-ring polyphenyl ether (MCS 210), fluorosilicone 
(XF 10301), and synthetic paraffin (XRM 109 F-3) from soft EHL 
model. 



Figure 10. — EHL lubrication coefficient as function of Bond number 
for four-ring polyphenyl ether (MCS 210), fluorosilicone (XF 10301), 
and synthetic paraffin (XRM 109 F-3) from soft EHL model. 
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Figure 11.— Angular position predicted by soft EHL model with 
C D = 0.52 and C EHL = 0.96 compared with data on fluorosilicone 
(MCS 210). Chamber pressure, 34.5 kPa; droplet radius, 1 mm; Bond 
number, 0.45; Weber number, l<We£j<3.2; Reynolds number, 
250<Re D <450. 


pressures. The data are the average C EHL over the three 
pressures at a given Bond number. The error bars reflect 
the uncertainty in the coefficient, the 68 percent 
probability that another set of measurements of velocity 
versus 9 of a given radius droplet of known properties for 
three chamber pressures would produce an average 
lubrication coefficient C EHL that would fall within the 
interval shown. Those data points without error bars 
were taken at one pressure and do not constitute an 
average. 

Figures 11 and 12 indicate the ability of the soft EHL 
model to predict measurements of angular positions as a 
function of velocity for a 1-mm-radius droplet of MCS 
210 at 34.5-kPa chamber pressure with C EHL of 0.96 and 
C D of 0.52. Figure 13 shows the predicted minimum film 
thickness as a function of 9 from the soft EHL model for 
the same droplet, chamber pressure, C D , and C EHL . 
Since h m[n was not measured directly, no data com- 
parison is shown in figure 13. 


Discussion 

The semiempirical physical models of a droplet of 
liquid supported on a gaseous film were corroborated by 
experimental measurements on three fluids in a rotating 
outer race in the presence of air at three chamber 
pressures. The unknown constants resulting from the soft 
EHL model (a drag coefficient C D and a lubrication 
coefficient C EHL ) were found to have plausible values 
(0.4 and 1.0, respectively) for those cases where surface 
tension forces dominated over gravity forces and air 

inertial forces. 

The experimental apparatus was suitable for testing the 
models as long as the droplets were not at excessively high 
angles. At high 9 the outer-race velocity was large enough 
to make the Weber number excessive, and the droplet 
shape deviated substantially from spherical. The film 
thickness became large because the radial force became 
smaller as 9 increased toward 90°, and the outer-race 
velocity was large enough to drag the droplet up to higher 
9 . The film Reynolds numbers here were too large to 
justify the neglect of the inertial terms in the film region. 
Experimentally the droplets were observed to oscillate in 
the y direction at angles above 60° and to stand off the 
race sufficiently to blow axially ( y direction) off the race 
and splatter against the viewing port glass. The nature of 
the outer-race geometry was such that it constituted a 
viscous pump, taking air in axially and driving it radially 
and tangentially out to the outer race and then 
tangentially and axially over the outer race back out to 
the chamber. The outer-race radius, r by =-6 . 55 mm, 
was necessary to keep the droplet from being carried out 
of the outer race axially. At the higher angles the film 
thickness became large enough to permit the droplet's 
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Figure 12.— Angular position as function of outer-race velocity 
predicted by soft EHL model with C D = 0.S2 and C EHL = 0.96 
compared with data on fluorosilicone (MCS 210). Chamber pressure, 
34.5 kPa; droplet radius, 1 mm; Bond number, 0.45; Weber number, 
1 <We i) <3.2; Reynolds number, 250 < Re/) <450. 



Figure 13.— Minimum Film thickness as function of angular position 
predicted by soft EHL model for C D = 0.52 and C EHL = 0.96 for 
fluorosilicone (MCS 210) chamber pressure, 34.5 kPa; droplet radius, 
1 mm; Bond number, 0.45; Weber number, l<We/)<3.2; Reynolds 
number, 250<Re D <500. 


escape when combined with the high air velocities driven 
by the large outer-race velocities. 

At the low end of the droplet angle spectrum the outer- 
race velocity was low as long as the droplet size was not 
excessive and the minimum thickness was small. This 
experimental domain was ideal for testing the theoretical 
models, the Weber numbers were small, and the film skin 
friction terms were becoming comparable to the air 
inertial drag terms so that differences in film lubrication 
models could have an appreciable effect on the results. 
When the air inertial drag dominated over the film skin 
friction, the tangential force balance was dominated by 


the balance of gravity forces and air drag, and the 
method of determining film thickness and the effect on 
the skin friction term in the tangential force balance 
became less significant. Practically it was difficult to keep 
the droplet at angles below 15° because the film thickness 
was becoming so small that the rms surface roughness 
combined with minor vibrations in the outer-race 
rotations caused asperity contact and the droplet wet the 
surface. 

Both models indicated that the contact area in the 
HL model and r}D x D y /4 in the EHL model) was 
independent of 6 . However, the EHL elliptical contact 
area was typically 20 percent larger than the HL contact 
area (i.e., D x D y /4r^=^\,2), The droplet was settling into 
the outer-race geometry more than was visually apparent 
for the low-Bond-number cases. 

No appreciable motion was observed in the oil droplets 
tested since the ratios of oil viscosity to air viscosity were 
typically 104. The neglect of these effects in the physical 
models appeared to be justified for these experimental 
conditions. 


Conclusions 

The dominant physical parameters governing the 
failure of lubricant droplets to wet nonconformal 
conjunctions because of a gaseous film between the 
droplet and the solid surface have been delineated by 
comparing two film lubrication models with experimental 
measurements. For sufficiently low Bond numbers 
PoixStd/o the experimental droplet angular position 0 on a 
rotating surface correlated well with two dimensionless 
groups: p a i r V 2 /poifdS and r) a i r V/p on r%g. Two constants 
were determined from experiment (a drag coefficient 
C D =z 0.4 and a lubrication coefficient C L - 3.0) when 
using a simple hydrodynamic model. A more complex 
soft elastohydrodynamic lubrication model correlated 6 
with the parameters P a wV 2 /p oi{ r d g, v aiT V/p oi[ r^g, r b>x /r d , 
and r by /r d for low Bond numbers with a drag coefficient 
C d ssQA and a lubrication coefficient Cehl — 1*0. 
Furthermore the soft EHL model predicted 0 for a wider 
range of Reynolds and Bond numbers if the model 
incorporated a drag-coefficient-versus-Reynolds-number 
relation and a lubrication-coefficient-versus-Bond- 
number relation. 

It has been suggested (Kingsbury, 1969) that, as the air 
pressure is reduced within the chamber, wetting of the 
outer race by the droplet occurs because the droplet- 
outer-race conjunction becomes starved of lubricant (i.e., 
air). However, lowering the air pressure only served to 
reduce the density of the air and hence the aerodynamic 
drag on the droplet in this configuration. The droplet 
remained fully flooded with lubricant (i.e., surrounded 
by air whose viscosity remained virtually unchanged). It 
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is true that as the pressure was reduced, the droplet sank 
to lower angles because of the reduction in aerodynamic 
drag. At the lower angles the component of gravity in the 
radial direction became greater and the film became so 
thin that asperity contact caused droplet wetting of the 
outer race. If, however, the velocity of the outer race was 
progressively increased as the chamber pressure was 
decreased, the droplet remained at a constant angular 
position and the film did not become so thin as to allow 
asperity contact and subsequent wetting of the outer race 
by the droplet. 
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